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ABSTRACT 


A  novel  airborne  GMTI  radar  approach  based  on  adaptive  multichannel  SAR 
processing  is  proposed  to  improve  the  surface  surveillance  sensor  performance  in 
detecting  slowly  moving  targets  on  the  ground.  The  key  sensor  attributes  are  broad 
bandwidth,  large  antenna  array  baseline  of  multiple  phase  centers,  and  adaptive  array 
processing  over  a  long  coherent  integration  time  interval.  This  sensor  architecture  may 
improve  the  GMTI  sensor  capabilities  in  several  aspects 

1)  Target  detection  at  very  low  minimum  detectable  velocity  (MDV) 

2)  Robust  adaptive  processing  to  cancel  strong  ground  clutter 

3)  High  sensitivity  to  detect  weak  targets 

4)  Flexible  array  requirement  including  sparser  arrays 

5)  Compatibility  with  SAR  imaging  applications. 


The  presented  sensor  architecture  is  suitable  for  deployment  on  a  fast,  stand-off  airborne 
surveillance  platform  against  weak  and  slow  targets  such  as  dismounted  troops, 
guerrillas,  terrorists,  and  illegal  vehicle  or  boat  traffic,  both  in  the  clear  and  under 
concealment. 
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1.  INTRODUCTION 


This  report  describes  a  novel  radar  and  processing  concept  that  addresses  wide-area,  standoff 
surveillance  and  detection  against  dismounted  troops,  guerrillas,  terrorists,  or  illegal  vehicle  or  boat 
traffic,  both  in  the  clear  and  under  concealment.  What  makes  this  problem  difficult  comes  from  a 
combination  of  factors. 


1 .  Targets  of  interest  listed  above  are  generally  small,  may  have  very  small  Doppler  signatures, 
and  may  be  concealed  in  trees  or  in  terrain  clutter. 

2.  Large  standoff  ranges  are  desirable  to  keep  the  surveillance  as  covert  and  as  safe  as  possible. 

3.  Areas  to  be  searched  can  be  very  large  and  large  area  surveillance  rates  are  needed. 

4.  To  catch  the  enemy  or  the  targets  in  action  or  in  exploitable  conditions,  persistent  coverage  is 
needed  which  requires  all  weather,  day,  and  night  sensing  capability. 


Current  microwave  surveillance  systems  can  provide  all  weather,  day  and  night,  rapid  and  wide  area 
coverage  surveillance  capabilities  against  vehicle  size,  and  stationary  or  moving  targets  in  clear  but  will 
be  limited  in  performance  against  the  small,  slow  and  concealed  targets  listed  above.  Synthetic  aperture 
radar  (SAR)  is  commonly  used  for  stationary  target  detection  and  space-time  adaptive  processing  (STAP) 
ground  moving  target  indicator  (GMTI)  methods  have  been  developed  for  moving  target  detection. 
Standard  SAR  use  a  single  phase-center  antenna,  wideband  waveforms,  and  long  coherent  integration 
times  to  achieve  high  resolution.  Moving  targets  typically  appear  smeared  in  SAR  images.  While 
algorithms  have  been  developed  [1-3]  to  refocus  target  returns  in  SAR  images,  a  single  phase  center  SAR 
is  inherently  limited  in  clutter  cancellation  due  to  lack  of  spatial  degrees  of  freedom  (DoF).  STAP  GMTI, 
on  the  other  hand,  uses  multiple  phase  centers,  which  provide  the  spatial  DoF  needed  for  clutter 
cancellation,  but  uses  narrowband  waveforms,  and  short  coherent  integration  times,  which  limit 
performance  against  very  small  and  slow  targets. 

This  report  describes  a  unified  multichannel  wideband,  long  CPI  (coherent  integration  interval) 
radar  concept,  USASTAR  (unified  synthetic  aperture  space  time  adaptive  radar),  for  fixed  wing  aircraft 
such  as  Global  Hawk,  Predator,  or  manned  aircraft,  that  combines  the  advantages  of  both  the  high 
resolution  SAR  and  multichannel  STAP  techniques  to  fill  in  the  technology  gap  mentioned. 

By  unifying  the  SAR  and  GMTI  properties,  this  radar  concept  offers  a  potential  capability  to  detect 
stationary  and  moving  ground  targets  at  arbitrary  speed  and  geometry.  The  USASTAR  technique  can  be 
applied  to  radar  operating  at  any  frequency  band  from  Ku  to  VHF,  though  in  this  report  only  X-band  and 
UHF  band  applications  are  presented.  While  the  X-band  application  can  provide  the  desired  higher 
resolution  for  target  discrimination  when  targets  are  in  the  clear,  the  UHF  application  provides 
penetration  through  some  concealment  such  as  foliage,  camouflage,  and  certain  kind  of  walls. 
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2.  THE  USASTAR  CONCEPT 


Detection  of  moving  targets  by  airborne  radar  is  generally  accomplished  by  multiple  phase 
centers  and  STAP  to  cancel  clutter  which  spreads  in  Doppler  due  to  the  aircraft  motion.  As  target  velocity 
becomes  smaller,  achieving  acceptable  performance  requires  very  large  antennas,  very  slow  or  stationary 
platforms,  or  much  higher  resolution  to  reduce  clutter.  An  approach  used  today  and  in  the  past  is  to  install 
a  pulse-Doppler  radar  on  a  stationary  (or  nearly  stationary)  platform.  The  new  concept  proposed  here  is 
based  on  a  faster  flying  aircraft,  providing  larger  standoff  range  and  larger  surveillance  rates,  with  a  very 
long  antenna  baseline,  and  higher  resolution. 

USASTAR  is  based  on  adaptive  clutter  cancellation  and  coherent  change  detection  using  multiple 
phase  centers  on  the  same  platform  (or  multiple  platforms  connected  with  a  high  precision  geotracking 
system  such  as  laser  geodimeter),  each  producing  a  SAR  image  or  a  high-resolution  range-Doppler  map 
data,  where  the  time  it  takes  for  the  platform  to  travel  the  distance  between  the  phase  centers  is  sufficient 
for  targets  to  move  a  number  of  cells  (USASTAR  can  be  used  even  when  targets  move  a  fraction  of  a  cell 
using  target  phase  change,  but  with  more  target  loss).  USASTAR  is  an  approach  to  GMT1  when  target 
Doppler  is  expected  to  be  very  low.  It  also  provides  simultaneous  capability  for  SAR  imaging  and 
detection  of  stationary  targets  and  GMTI  detection  of  moving  targets.  The  concept  is  illustrated  in 
Figure  1 . 


Figure  1.  USASTAR  concept. 


It  is  expected  that  USASTAR  will  provide  improved  clutter  cancellation  and  sensitivity,  as 
compared  to  existing  SAR  change  detection  or  GMTI  systems,  in  detecting  small  and  slow  movers  such 
as  dismounts.  USASTAR  can  achieve  this  improvement  over  conventional  SAR  pass-to-pass  change 
detection  by  providing  better  registered  images  from  multiple  phase  centers  on  the  same  platform  and  by 
reducing  the  time  between  the  radar  looks  (depending  on  the  length  of  the  baseline  and  the  speed  of  the 
aircraft),  thereby  reducing  the  uncorrelated  clutter  returns  between  the  various  looks.  Because  USASTAR 
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produces  multiple  channel  high-resolution  SAR  or  range-Doppler  map  data,  more  clutter  sample  data  are 
available  for  local  adaptation  and  training  of  adaptive  array  processing  algorithms  used  in  GMTI 
processing.  Thus,  the  USASTAR  approach  enables  more  accurate  in-situ  data-driven  array  calibration, 
better  mitigation  of  array-airframe  interactions,  and  higher  cancellation  of  strong,  inhomogeneous  clutter 
in  realistic  field  environments,  compared  to  conventional  GMTI  systems. 

In  the  following  sections  two  alternative  formulations  of  the  USASTAR  concept  and  predicted  results 
are  presented.  Also  presented  is  a  sample  result  from  application  of  the  method  to  Lockheed  Martin 
Tuxedo  X-band  SAR  data  [4]. 
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3.  USASTAR  FORMULATION  AND  PERFORMANCE  PREDICTIONS 


Theoretical  tools  have  been  developed  to  help  define  optimal  USASTAR  system  architecture  and 
parameters.  The  key  USASTAR  system  attributes  are  broad  bandwidth  large  antenna  array  baseline,  and 
long  coherent  integration  time  (CPI).  Both  bandwidth  and  array  structure  are  important  system  design 
parameters  that  need  to  be  carefully  selected  to  trade  off  the  intended  system  performance  goal  in 
operational  utility  and  functionality  against  the  system  complexity  and  cost.  These  two  parameters 
significantly  affect  radar  system  performance  in  suppressing  radio  interference,  jamming,  and  clutter. 
Furthermore,  long  CPI  enhances  radar  Doppler  (or  cross-range)  resolution  that  is  essential  to  increase 
radar  target  detection  sensitivity  and  clutter  suppression  performance.  For  GMTI  radar,  the  ability  to 
cancel  strong  clutter  is  a  prerequisite  to  reliable  detection  of  slowly  moving  targets.  In  radar  design,  the 
trade-off  analysis  of  adaptive  array,  signal  processing,  and  performance  typically  relies  on  numerical 
simulation.  While  computer  simulations  may  be  useful  in  some  circumstance,  this  approach  may  demand 
substantial  computational  resource  to  simulate  a  system  characterized  by  a  large  bandwidth  and  a  large 
array  of  many  DoFs.  In  the  latter  situation,  it  will  be  very  desirable  to  have  analytic  tools  that  can 
accurately  and  efficiently  predict  the  system  performance  in  clutter  cancellation.  For  this  purpose,  this 
section  will  present  a  set  of  analytic  formula  that  may  be  used  to  compute  the  radar  clutter  cancellation 
performance  in  terms  of  a  given  radar  waveform  in  its  signal  frequency  composition  and  bandwidth  and 
of  an  array  structure  in  its  element  placement  and  aperture. 

3.1  SNR  LOSS  IN  LONG-CPI  PROCESSING 

A  set  of  closed-form  expressions  have  been  derived  that  specify  the  theoretical  lower  bound  of 
the  signal-to-noise  ratio  (SNR)  loss  incurred  in  the  clutter  cancellation  via  adaptive  array  processing.  The 
radar  parameters  and  assumptions  used  in  the  theoretical  development  are  defined  as  follows. 

Here,  the  velocity  of  the  airborne  radar  aircraft  is  assumed  to  be  va  moving  along  a  straight  line 
parallel  to  a  linear  radar-receiving  array.  The  linear  array  consists  of  N  array  sensor  elements  whose 
phase  centers  are  located  at  D„,  n  =  0,  1,  2,  ...  N  -  1  as  measured  from  a  reference  location  such  as  the 
phase  center  of  a  separate  radar  transmitter  antenna.  The  placement  of  the  N  sensor  elements  can  be 
arbitrarily  distributed  within  the  linear  array  aperture  D  that  may  be  sparse  or  only  partially  populated.  In 
order  to  minimize  the  SNR  loss  in  adaptive  clutter  cancellation,  radar  target  signal  integration  will  be 
carried  out  over  a  long  CPI.  In  this  limit,  the  CPI  will  be  longer  than  the  time  of  radar  flight  over  the 
entire  array  aperture  length  D,  CPI  »  D/v„.  As  target  migration  is  expected  over  a  long  CPI  due  to 
range-walk  of  a  moving  target  and  its  signal  smearing  cross  multiple  Doppler  filters,  proper  target  motion 
compensation  will  be  assumed  to  avoid  any  target  signal  loss  incurred  in  integration.  The  radar  waveform 
is  characterized  by  a  spectrum  pg(v)  that  spans  a  frequency  band  [-B/2.  B/2]  center  at  fo.  The  signal 
spectrum  pn(v)  may  be  arbitrary  and  sparse  in  its  spectral  occupation  within  the  fractional  bandwidth  Bf  = 
B/f).  However,  the  spectrum  is  known  such  that  perfect  range  compression  or  matched  filtering  of  the 
pulse  signal  is  feasible. 

Given  these  parameters,  the  SNR  loss  at  the  target  radial  velocity  v  incurred  by  an  optimal  array 
processor  in  nulling  clutter,  that  is  characterized  by  the  wind-blown  internal  clutter  motion  (ICM) 
spectrum  Pk-m(v),  can  be  shown  to  be 
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SNR  loss  =  tHRrJy  (v)t  ; 

Rcn{v)  =  (cnr)02  [_Picm{u-v)z{u)zH  {u)du  +  \N  ; 


Target  vector  [t]  (  =-y=  ,  clutter  vector  [c(v)]  = 


h{v) 


Here,  the  t  is  the  target  steering  vector  and  Rcn(v)  is  the  clutter-plus-noise  covariance  matrix  that 
can  be  expressed  as  the  convolution  of  the  wind-blown  clutter  ICM  spectrum  Picm(v)  and  the  outer 
product  of  c(v).  The  latter  quantity  is  actually  the  clutter-plus-noise  covariance  in  the  absence  of  any 
clutter  spread  due  to  internal  clutter  motion  and  is  a  function  of  <p(v). 


$,(v)  =  exP (-M/O  {  pB{v)cxpi-j(odTnBf^\dv  ,  ^  =  £|4,(v)| 


2K  D  B 

0)Q  —  2n ,  k0  —  — —  ,  cod  —  2k0v  ,  Tn  —  —  ,  B ^  —  ——  . 

K  /o 

The  function  </>(v)  is  related  to  the  Fourier  transform  of  the  waveform  spectrum  ps(vj  and  its 
functional  forms  for  two  special  cases,  the  constant  and  the  Gaussian  spectra,  are  given  as  follows 


t  =exp(-/^/r„) 


sin  (codr„  Bf/2) 
a)dTn  Bf/2 


for  constant  spectrum  pR  ( v)  = 


=  exp  (- j 6) dTn ) exp  -  -  dT"Bf 

Qtt 


2  (  Av2 

for  Gaussian  spectrum  ps(v)  =  —*=  exp - — 

v  7tB  l  B~ 


The  clutter  ICM  spectrum  Picm(v)  is  composed  of  two  components,  the  DC  part  S(v)  and  the  AC  spectrum 
Pac(v)  due  to  motion  modulation: 

V  1 

Clutter  ICM  spectrum p1CM  (v)  = - -J(v)  +  - — - -pAC(y)  ; 


Clutter  correlatoin  plCM  (t)  =  p/CM  (v)exp(j2k0vr)dv  . 


One  popular  model  is  to  employ  an  exponential  spectrum  or  its  corresponding  temporal  correlation 
function  piCM(P  to  characterize  the  clutter  ICM 
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Exponential  ICM  spectrum:  pAC  (v) 


2fl 


■exp 


f  v1^ 


A 


v  >  >•  y 


Clutter  correlatoin:  pICM  (r)  = 


-  +  - 


1 


1 


r  +  \  r  +  U  +  4k2p2T2 


DC/AC  power  ratio  r: 

101oglo  r  =  63.2  -12.1  logl0( Frequency,  MHz) 

- 1 5.5 logl0 (Wind  speed,  mph) 

Clutter  spread:  J3V  =  0.1048(log10(wind  speed,  mph) +  0.4147]  m/s  . 


(3.5) 


In  Equation  (3.5),  both  functional  parameters,  the  DC/AC  power  ratio  r,  and  the  root  mean  square 
(rms)  spectral  spread  /?,,,  are  dependent  on  the  surface  wind  speed.  Unfortunately,  the  exponential  clutter 
ICM  model  does  not  render  Equation  (3.1)  into  a  mathematically  tractable  form.  Thus,  for  the  reason  of 
obtaining  a  mathematical  formula  convenient  for  computation,  an  equivalent  Gaussian  ICM  model  is 
sometime  preferred.  This  model  is  defined  by 


Gaussian  ICM  spectrum:  pAC  (v)  =  .  .  exp  - 

V2*#2  l 

Clutter  correlation:  pAC  (t)  =  exp(-2 k2p2r2 


,.2  A 


2  PI) 
)  ■ 


(3.6) 


The  clutter  spread  parameter  /?„  for  the  Gaussian  ICM  spectrum  typically  must  assume  a  much 
larger  value  than  the  corresponding  /?,,  used  in  the  exponential  model  in  order  to  result  in  comparable  SNR 
loss  results  as  predicted  by  Equation  (3.1).  By  assuming  both  a  Gaussian  clutter  ICM  and  a  Gaussian 
waveform  spectra,  a  closed  form  expression  for  the  clutter-plus-noise  covariance  function  Rcn(v)  of 
Equation  (3.1)  can  be  obtained  for  this  special  case  as 


r  + 1 


c(v)cw(v)]k 


r  + 1 


G . = 


1  + 


([c(v)cw(v)]mn/7,r(r„-Tj)' 
0.5(Bfkopv)2(T2m  +  T2nj\  . 


+  S.„ 


Also,  the  following  approximation  of  Rcn(v)  is  useful  for  narrowband  waveforms 


(3.7) 


R cn  (v)  =  (CNR)<f>2  £ p,CM  ( u-v)c(u)cH  (u)du  +  1N 

=  (CNR)<f>2c(v)cH  {v)  O  R/ow  +I/V  . 
for  k0pv Bf  max (rn)«l  ; 

[R/CA/  =  PlCM  (TI1  ~^m)  ■ 
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In  the  absence  of  any  wind-blown  clutter  modulation.  Equation  (3.1)  may  be  considerably  simplified  as 


SNR  loss  =  l  — 


(< CNR)(fr 


\  +  (CNR)(/)2 


t"c 


(3.9) 


In  this  situation,  the  SNR  loss  may  be  obtained  from  Equations  (3.2)  and  (3.3)  by  substituting  the 
following  inner  product  of  t  and  c  into  Equation  (3.9): 


7H-  1  v*  /  ,  xSin(^5/„) 

t  c  =  --£cxp{-j2irfy„)  — 

■/V  n= 0 


for  constant  spectrum  pB  (v)  =  —  ,  ^  ,  C  is  the  light  velocity. 

B  Cva 

1  N-l  f  02.2") 

twc  =  —  Z  exP  (  -./>'2^/o4  )  exp  — -f- 

2V  «=  0  V  4  / 

2  f  4t>2  ^ 

for  Gaussian  spectrum  pB(v)-  r—  exp 

yj/rB 


(3.10) 


V  B2  J 


These  equations  also  may  be  modified  to  accommodate  a  waveform  spectrum  that  may  span  several  non¬ 
overlapping  subbands  within  the  system  bandwidth  B.  For  example,  if  there  are  M  separate  subbands  b,„ 
individually  centered  at  f„„  m  -  0,  1,  ...  M  - 1  within  B,  each  with  a  constant  power  spectral  density  1/BM, 
Bm  ~  b0  +  b /  ...  +  bu-i,  then  Equation  (3.10)  becomes 


v^1  1  ,  ..  -  .  sin  (nb  t  ) 

m= 0  n=o 


M- 1 

Pb{v)=1LP’M  »  PAv)=~ir  .  w  =  0,l,  ...  M-l  and 

m=0 


1 


5, 


(3.11) 


ay 


A/ -I 


«'®[/.-0-»../.+0.5*.]e[/„-#.5«,/0+0.5B].JK  =  2;».  • 


m=0 


3.2  COMPUTATIONAL  EXAMPLES 

Theoretical  SNR  loss  due  to  the  clutter  residue  as  the  output  of  an  imperfect  clutter  nulling 
process  are  plotted  in  Figures  2  through  4  for  a  few  sample  cases.  These  figures  illustrate  the  SNR  loss 
dependence  on  array  structure,  bandwidth,  radar  speed,  clutter  strength,  and  wind  speed.  Generally 
speaking,  a  low  SNR  loss  indicates  good  clutter  cancellation  and  better  adaptive  array  processing 
performance  that  results  in  a  lower  minimal  discernible  detectable  velocity  of  a  moving  target.  By 
inspecting  these  figures,  it  is  clear  that  the  best  results  are  achieved  by  using  an  array  with  a  long  aperture 
baseline  and  a  large  bandwidth  to  suppress  the  periodic  loss  due  to  the  grating  lobe  effects  when  the  array 
is  only  populated  by  sparse  elements  distributed  over  a  large  aperture. 

Figure  2  compares  the  SNR  loss  of  an  array  with  two  separate  subarray  apertures  operated  with 
either  the  narrowband  mode  shown  in  red  in  Figure  2(a),  or  the  broadband  mode,  shown  in  green  in 
Figure  2(b),  with  a  bandwidth  of  200  MHz  centered  at  300  MHz.  Two  array  structures  are  illustrated. 
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The  first  is  a  filled  array  of  101  elements  spanning  a  total  array  length  of  50  m.  The  second  is  a  sparse 
array  consisting  of  two  subarrays  separated  by  40  m.  Each  subarray  has  1 1  antenna  elements  in  half¬ 
wavelength  spacing  over  a  span  of  5  m.  For  reference,  the  SNR  loss  curves  for  a  filled  array  are  also 
plotted;  one  is  for  the  case  of  a  stationary  radar  (in  black)  whereas  the  other  an  airborne  radar  moving  at 
the  speed  of  50  m/s  (in  blue).  In  these  charts,  the  solid  curves  are  modeled  by  assuming  an  Gaussian 
clutter  ICM  spectrum  for  a  very  windy  condition  at  the  speed  of  20  mph  such  that  the  DC-to-AC  ratio  r  is 
11.1  dB  and  the  rms  clutter  spread  /?,.  is  0.4  m/s.  The  ideal  cases  in  the  absence  of  wind  are  plotted  in 
dashed  curves. 

Figure  3  compares  the  SNR  loss  of  two  sample  cases.  The  first  case  (in  red)  is  generated  for  a 
short  array  aperture  of  4  m  in  length  and  is  operated  over  a  20  MHz  bandwdith  between  380  and  400 
MHz.  While  the  second  case  (in  green)  assumes  two  subarray  apertures,  each  of  2  m  in  length,  separated 
by  500  m  but  the  radar  band  spans  from  200  to  400  MHz.  The  radar  speeds  for  both  cases  are  taken  to  be 
50  m/s.  Due  to  the  different  bandwidth  and  consequently  different  range  resolution,  the  element  clutter  to 
noise  ratio  (CNR)  of  50  dB  for  the  first  case  is  10  dB  stronger  than  that  of  the  second  case.  Also,  the 
wind-blown  clutter  ICM  spectrum  is  taken  as  Gaussian.  At  the  surface  wind  speeds  of  5  mph  and  20 
mph,  respectively,  the  ICM  DC-to-AC  ratios  are  20  dB  and  11.1  dB,  respectively,  and  the  rms  clutter 
spread  are  0.2  m/s  and  0.4  m/s. 

Figure  4  compares  the  difference  in  SNR  loss  for  two  clutter  levels.  In  this  exercise,  the  radar 
moves  at  a  much  faster  speed  of  200  m/s.  The  array  consists  of  two  small  subarrays  of  two  wavelengths, 
or  five  elements  separated  by  500  m.  The  radar  frequency  band  spans  from  200  to  400  MHz.  The  same 
wind-blown  clutter  ICM  model  is  used  as  that  for  Figure  3.  As  clearly  shown  by  this  figure,  the  SNR  loss 
in  the  presence  of  a  lower  clutter  level  at  the  element  CNR  of  20  dB  are  clearly  smaller  than  the 
corresponding  cases  at  the  element  CNR  of  40  dB. 

3.3  SAR-BASED  ANALYSIS  AND  COMPUTATIONAL  EXAMPLES 

In  this  section,  an  alternative  formulation  of  the  USASTAR  processing  concept  and  predicted 
performance,  as  measured  by  signal-to-interference  noise  ratio  (SINR)  loss  factor,  will  be  presented.  In 
this  formulation,  the  long  CPI  wideband  Doppler  processing  will  be  treated  as  SAR  processing,  and  the 
USASTAR  concept  will  be  treated  as  adaptive  processing  of  multichannel  SAR  images  [5].  Effects  due 
to  clutter  internal  motion  are  not  included  here,  but  could  be  added  along  the  same  lines  as  in  the  previous 
analysis. 

Figures  5  through  14  show  the  properties  associated  with  SAR  processing  prior  to  spatial  adaptive 
processing,  theoretical  formulas  for  SAR  and  USASTAR  steering  vectors  and  predicted  SINR  loss 
performance,  under  simplifying  assumptions  shown,  compared  to  conventional  GMT1.  Both  X-band  and 
UHF  examples  are  given. 
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Radar  speed:  50  m/s 
Element  CNR:  40  dB 


Target  Range  Rate  (m/s) 


I - 1 - 1 

0.01  0.1  1 

Target  Range  Rate  /  Radar  speed 

(a) 


I - 1 - 1 

0.01  0.1  1 

Target  Range  Rate  /  Radar  speed 

(b) 


Center  frequency:  300  MHz 
Bandwidth 

Narrowband  (black  and  blue) 

Linear  array  phase  centers 

1)  101  at  DO  :  0.5  : 50  ]  X 
(stationary  platform) 

2)  101  at  [0  :  0.5  :  50  ]  X 

3)  22  at  [0  :  0.5  :  5  ] 
and  [45  :  0.5  :  50]  X 

Clutter  ICM  model 
Dashed  line:  no  wind 
Solid  line 

20  mph,  Gaussian  spectrum 
DC/AC:  11.1  dB 
RMS  spread:  0.4  m/s 


Radar  speed:  50  m/s 
Element  CNR:  40  dB 

Center  frequency:  300  MHz 
Bandwidth 

Narrow  band  (black  and  blue) 
200  MHz  (green) 

Linear  array  phase  centers: 

1)  101  at  [0  : 0.5  :  50  ]  X 
(stationary  platform) 

2)  101  at  [0  :  0.5  :  50  ]  X 

4)  22  at  [0  :  0.5  :  5  ] 

and  [45  :  0.5  :  50]  X 

Clutter  ICM  model 
Dashed  line:  no  wind 
Solid  line 

20  mph,  Gaussian  spectrum 
DC/AC:  11.1  dB 
RMS  spread:  0.4  m/s 


Figure  2.  SNR  array  loss  compared  with  two  separate  svbarray  apertures  operated  with  (a)  narrowband  mode,  or 
(b)  broadband  mode. 
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SNR  Processing  Loss  (dB)  SNR  Processing  Loss  (dB) 


Target  Range  Rate  (m) 


Platform  speed:  50  m/s 
Long  CPI  (»D/V~10s) 

Clutter  ICM  model: 
Gaussian  spectrum 
Wind  speed:  0,  5,  20  mph 

Array  structure 

*  Single,  filled  array: 

-One 4m  aperture 

-  Band:  380-400  MHz 
-CNR:  50 dB 

•  Two  sparse  arrays: 

-  Two  2  m  apertures, 
500  m  separation 

-  Band:  200-400  MHz 
-CNR:  40  dB 


Figure  3.  SNR  array  loss  for  two  sample  cases. 


Radar  speed:  200  m/s 
Frequency  band:  200-400  MHz 

Linear  array: 

2  small  filled  apertures, 
each  has  5  elements 
Element  spacing:  0.5  X 
Aperture  separation:  500  m 

Array  element  CNR: 

20  and  40  dB 
Clutter  ICM  spectrum: 
Gaussian  spectrum 
Wind: 

Solid  line:  No  wind 
Dotted  line:  5  mph 
Dashed  line:  20  mph 


Figure  4.  SNR  array >  loss  for  two  clutter  levels. 
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Space  (phase  centers) 


Data  Cube  Vectorized  Linear  transformation 

Data  Cube  to  SAR  basis 


Multichannel  SAR 
image 


>s1 

♦  s2  =  M2z2 
>  s3  =  M3Z3 

-►  s  4  =  M4z4 


Cross-range 


Stagger  data  in 
slow  time  to  meet 
DPCA  condition 


i  SAR  processor 
knows  locations  of 
phase  centers 

Each  SAR  resolution  cell 
is  automatically  “phased 
up”  for  stationary  clutter 
at  that  location 

Figure  5.  USA  STAR  SAR  preprocessing  concept. 


•  SAR  is  the  wideband,  long-CPI  generalization  of  the 
Doppler  processor  in  ordinary  post-Doppler  STAP 

•  Linear  transform  of  input  data  cube 

-  Annihilates  the  exoclutter  subspace 

-  Invertible  transformation  of  the  endoclutter  subspace 

*  “Freezes”  clutter  into  SAR  resolution  bins 

-  Clutter  in  one  bin  is  well-decorrelated  from  other  bins 
multiple  resolution  cells  away. 

-  Stationary  targets  have  trivial  steering  vectors 

-  Stationary  clutter  has  trivial  covariance 

*  Moving  targets  smear  over  multiple  resolution  cells 

Figure  6.  Properties  of  SAR  preprocessing. 
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Slant  Range 


Image  Domain 


"(Xs.Ys) 


Cross  Range 


Spatial  Frequency  Domain 


Li 


Cross  Range  Spatial  Frequency 


Figure  7.  Theoretical  formula  for  SAR  steering  vector. 


Other  Elements  Transmit  Phase  Center  Receive  Element  of  Interest 
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Cross-track  Ground  Velocity  {m/s) 


*  Large  clutter-to-noise  ratio 

*  Elements  are  mutually  calibrated 

*  No  internal  clutter  motion,  crab,  unmeasured  aircraft 
motion  and  vibration 

*  No  jammers  and  other  interference 

*  Isotropic  element  patterns 

*  Optimal  AMF  processing  with  perfect  knowledge  of  steering 
vectors  and  clutter  covariance 

Figure  9.  Simplifying  assumptions  for  theoretical  SINR  loss  performance  predictions. 
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Figure  10.  Theoretical  formulas  for  USASTAR  SINR  loss. 


-15  -10  -5  0  5  10  15 

dB 

Along-track  Ground  Velocity  (m/s) 


Carrier  =  1  GHz 
Bandwidth  =  10  MHz 
CPI  =  SO  ms 
Range  =  10  km 
Airspeed  =  180  m/s 
Phase  center 
locations  =  [0,1]  m 


Figure  1 1.  SINR  loss  estimate  for  conventional  GMTI. 
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Cross-track  Ground  Velocity  (m/s) 


•Target  at  10  km  slant  range 

•Loss  computed  using  optimal  processor  with  full 
knowledge  of  clutter  covariance 


|  *4  MHz  bandwidth 
;  -435  MHz  carrier 
5  *3  second  CPI 
j;  -130  m/s  platform  speed 
j  *3.74  m  linear  array 
» *12  elements 


s.  -30  MHz  bandwidth  i  -200  MHz  bandwidth 

■  -435  MHz  carrier  ;  -350  MHz  carrier 

l 

:  *3  second  CPI  •  -30  second  CPI 

|  t 

f  >50  m/s  platform  speed  j.  -180  m/s  platform  speed 

|  -10  m  linear  array  |  -30  m  sparse  array 

•30  elements  :  -Elements  at  [0,2,8,30]  m 


Figure  13.  Estimated  SINR  loss  for  three  different  sets  of  radar  parameters  at  UHF. 
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•Target  at  10  km  slant  range 

•Loss  computed  using  optimal  processor  with  full 
knowledge  of  clutter  covariance 


•200  MHz  bandwidth 
•350  MHz  carrier 
•30  second  CPI 
■•180  m/s  platform  speed 
•30  m  sparse  array 
•Elements  at  [0,2,8,30]  m 


•200  MHz  bandwidth 
•350  MHz  carrier 
•30  second  CPI 
•180  m/s  platform  speed 

! <30  m  sparse  array  +  tow 

[ 

•Elements  at  [0,2,8,30,300]  m 


Figure  14.  USASTAR  SINR  loss  predictions  showing  the  effect  of  extending  the  antenna  baseline  with  a  towed 
aperture  300  m  behind  the  first  antenna  element,  at  UHF. 


4.  HIGH-RESOLUTION  ADAPTIVE  ARRAY  PROCESSING  APPROACHES 


Airborne  GMTI  radar  operation  confronts  many  technological  challenges.  Most  of  all,  the  radar 
signal  processing  must  cancel  strong  ground  clutter  interference  with  a  large  Doppler  bandwidth  in  order 
to  detect  slowly  moving  ground  targets.  Typically,  such  a  task  is  accomplished  by  utilizing  a  receiving 
array  of  many  phase  centers  and  multichannel  receivers  that  are  followed  by  adaptive  array  processing. 
To  achieve  good  clutter  cancellation  performance  in  a  realistic  field  operational  environment,  the 
processing  algorithms  must  overcome  many  practical  difficulties.  For  example,  different  channels  of  the 
receiver  system  may  not  be  well  matched.  The  array  response  may  be  distorted  as  a  result  of  near-field 
array-airframe  interactions.  Furthermore,  the  clutter  reflection  from  the  ground  is  typically  not  only 
strong  but  also  heterogeneous  due  to  various  terrain  types  under  radar  illumination  and  various  discrete 
and  diffuse  surface  scatterers. 

Mathematically,  good  adaptive  array  processing  performance  requires  accurate  estimation  of  the 
clutter  covariance  matrix  R  and  the  actual  array  response  v  to  a  target.  In  practice,  R  is  substituted  by  its 
sample  covariance  matrix  computed  or  trained  with  sufficient  clutter  samples  that  should  be  many 
multiples  of  the  dimensionality  of  the  data  vector  to  avoid  excessive  estimation  loss.  These  training 
samples  must  be  statistically  identical  to  the  true  clutter  to  be  canceled  that  is  imbedded  in  the  test  data 
vector  2,  a  condition  very  difficult  to  satisfy,  especially  in  a  heterogeneous,  nonstationaiy  ground  clutter 
environment.  The  estimation  of  R  also  may  be  plagued  by  the  presence  of  other  interfering  targets  or 
target  like  signals  in  the  clutter  training  ensemble  that  cause  undesirable  mutual  target  interference. 
Furthermore,  the  phase  response  to  a  target  across  the  linear  receiving  array  will  be  distorted  due  to  the 
strong  near-field  multipath  interactions  between  the  array  and  the  airframe  or  other  anomalies,  including 
mismatched  frequency  response  among  different  radar  channels. 

Nonlinear  array  phase  distortion  is  a  serious  concern.  The  adaptive  array  processor  is  normally 
tuned  to  the  desired  target  signal  whose  array  response  must  be  precisely  known.  To  achieve  both  high 
signal  gain  and  clutter  cancellation,  one  must  calibrate  the  array  response  to  accurately  predict  the  array 
signal  attributed  to  a  target. 

In  this  section,  several  adaptive  algorithms  will  be  described  to  achieve  superb  clutter 
cancellation  and  GMTI  detection  performance  in  either  radar  range-Doppler  or  SAR  image  domain.  One 
algorithmic  approach  utilizes  an  adaptive  matched  filter  (AMF)  that  is  matched  to  the  array  response  to 
the  expected  target  as  derived  from  local  clutter  eigenvector.  This  algorithm  employs  a  “local”  array 
calibration  process  to  extract  moving  target  signal  in  the  endo-clutter  region  but  is  unable  to  derive  the 
array  target  response  in  the  exo-clutter  region  without  clutter  data.  A  similar  but  different  approach 
implements  the  single-cell  AMF  into  two  steps.  Initially,  a  “global”  array  calibration  matrix  used  for 
correcting  array  distortion  is  derived  from  radar  range-Doppler  clutter  data  via  a  least-square  procedure. 
Next,  the  calibrated  target  steering  vector  is  incorporated  in  the  AMF  to  achieve  high  array  SINR 
improvement  in  both  endo-  and  exo-clutter  regions.  A  different  category  of  more  complicated  multicell 
algorithms  capable  of  better  clutter  suppression  also  will  be  defined.  One  type  is  a  generalized  matched 
filtering  process  in  the  two-dimensional  space-fast-time  domain  that  processes  additional  fast-time  range- 
delay  data  samples  in  each  array  element  channel.  Another  technique  of  a  similar  approach  is  a  two- 
dimensional,  multiple-pixel  clutter  whitening  process  incorporating  a  noise-canceling  scheme. 

A  common  starting  point  of  all  these  algorithms  is  a  high-resolution  multichannel  data  set.  Given 
a  stream  of  range-pulse  data  from  each  receiving  array  channel,  a  range-Doppler  map  or  SAR  image  is 
initially  formed  by  integration  over  a  long  CPI.  Later,  several  or  all  of  the  data  maps  or  images  from 
multiple  array  channels  are  to  be  adaptively  combined  to  achieve  the  desired  clutter  cancellation  and 


17 


target  detection  results.  The  adaptive  array  processing  techniques  to  be  explored  here  may  be  viewed  as 
different  element-space,  post-Doppler  adaptive  array  beamforming  processes.  In  actual  implementations 
of  Doppler  filtering  or  SAR  imaging,  coherent  integration  is  performed  over  a  relatively  long  CPI.  The 
long  CPI  results  in  improved  radar  Doppler  or  cross-range  resolution.  Thus,  more  clutter  training  samples 
will  be  available  from  the  neighborhood  of  a  test  target  cell  to  better  estimate  clutter  covariance.  In 
practice,  the  number  of  local  clutter  training  samples  may  approach  100  times  the  data  vector  dimension. 
The  large  clutter  sample  size  facilitates  a  very  useful  data  target  excising,  which  is  a  selective  training 
procedure  found  to  be  very  effective  in  mitigating  the  mutual  target  interference.  High-resolution  range- 
Doppler  or  SAR  image  data  also  provide  sufficient  clutter  samples  to  support  local  array  calibration  that 
may  correct  array  response  distortion  in  fine  elevation  and  azimuth  angle  increments.  The  opportunity  for 
fine-grain  local  calibration  is  one  important  element  that  results  in  good  adaptive  array  processing 
performance. 

There  is  no  fundamental  distinction  between  radar  range-Doppler  and  SAR  image  processing 
except  for  the  motion  compensation  applied  to  radar  data.  The  SAR  imaging  is  a  process  that  maps 
ground  scatterers  onto  the  radar  range  and  cross  range  data  domain  by  applying  precise  aircraft  motion 
compensation  via  a  suitable  mapping  technique  such  as  the  backprojection  image  formation.  On  the  other 
hand,  fast  Fourier  transform  (FFT)-based  Doppler  filtering  is  commonly  applied  in  range-Doppler 
processing  to  radar  range-gated  data  in  order  to  transform  the  ground  onto  a  range-Doppler  map.  Because 
the  integration  is  carried  out  over  a  long  CPI,  appropriate  motion  compensation  also  may  be  necessary  in 
range-Doppler  processing.  The  compensation  could  be  only  approximately  and  reasonably  accurate  but 
not  necessarily  exact;  otherwise,  range-Doppler  processing  essentially  becomes  synonymous  with  SAR 
image  processing. 

4.1  SINGLE-CELL  ADAPTIVE  MATCHED  FILTERING  WITH  “LOCAL”  ARRAY 

CALIBRATION 

The  adaptive  processing  algorithms,  emphasized  in  this  report,  generally  belong  to  a  category 
known  as  element-space  post-Doppler  adaptive  beamforming.  In  particular,  the  Doppler  filtering  step 
prior  to  array  beamforming  is  applied  to  a  full  range  of  radar  pulse  data.  The  coherent  integration  is 
carried  out  over  a  relatively  long  CPI  to  increase  target  signal  gain  and  better  prepare  high-resolution 
radar  data  for  further  adaptive  processing.  The  long-CPI  Doppler  filtering  necessitates  several  essential 
procedures  in  the  signal  processing  chain  to  ensure  good  clutter  suppression  and  GMTI  target  detection. 
These  procedures  include  pulse  data  phase  correction  for  motion  compensation,  local  data-selective 
training  for  sample  clutter  covariance  estimation,  and  array  calibration.  Due  to  the  radar  aircraft  motion, 
the  Doppler  shift  of  a  scatterer  changes  appreciably  over  a  long  CPI.  To  correct  the  Doppler  drift, 
quadratic  phase  compensation  as  a  primitive  form  of  SAR  processing  is  applied  to  the  pulse  train  prior  to 
Doppler  filtering.  The  selective  training  is  a  procedure  that  iteratively  runs  the  same  adaptive  algorithms 
more  than  twice.  During  the  initial  iteration,  a  low  detection  threshold  is  set  to  select  all  potential  moving 
targets  in  the  test  data  output.  The  initial  detection  result  is  used  to  selectively  delete  all  potential  moving 
targets  from  the  clutter  covariance  training  data  in  subsequent  runs. 

The  local  array  calibration  procedure  may  be  integrated  as  part  of  the  post  Doppler  adaptive  array 
processing.  Given  an  array  data  vector  z  of  dimension  A  to  be  tested  for  the  presence  of  a  target  signal  v 
at  the  range  r  and  the  Doppler  /d,  one  initially  estimates  the  local  clutter  covariance  matrix  R  and  its 
principal  eigenvector  u,  then  the  AMF  output  power  SAMF  is  computed  by  using 
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In  Equation  (4.1),  the  notation  denotes  the  radar  wavelength  as  2,  the  array  element  spacing  as  d, 
the  depression  angle  as  0,  that  varies  with  the  target  range  r  and  the  aircraft  altitude  H  as  s\n~]  (H/r).  The 
azimuth  angle  of  radar  line-of-sight  (LOS)  directed  to  the  range-Doppler  cell  relative  to  the  array 
broadside  may  be  calculated  from  the  target  Doppler^,  the  aircraft  speed  va,  and  X  as  shown.  The  sample 
clutter  covariance  matrix  R  is  computed  from  sufficient  sample  data  vectors  x  selected  from  the 
neighborhood  of  the  range-Doppler  cell  z  under  test.  The  principal  eigenvector  u  of  R  captures  the 
distorted  array  response  specific  to  the  direction  as  defined  by  the  LOS  depression  angles  #and  0d.  The 
target  steering  vector  v  is  modeled  such  that  its  «th  vector  component  is  the  product  of  the  corresponding 
components  of  an  ideal  linear  array  response  vector  and  the  principal  clutter  eigenvector  u.  For  each 
beam  position  </>h  one  set  of  AMF  output  is  computed  for  all  range-Doppler  cells. 

It  is  intuitively  clear  that  the  post-Doppler  processing  procedures  can  be  readily  extended  to 
multichannel  SAR  processing  by  replacing  the  range-Doppler  map  for  each  element  channel  with  its 
equivalent  SAR  image  in  the  slant-range  and  cross-range  coordinates.  The  SAR  processor  accounts  for 
all  range  and  Doppler  migration  over  the  entire  CPI  due  to  along-track  aircraft  motion,  thereby  achieving 
higher  coherent  signal  gain  and  better  cross-range  image  resolution  from  along-track  integration. 
Usually,  the  radar  pulse  returns  for  each  channel  are  independently  processed  via  fast  backprojection  into 
a  complex,  co-registered  SAR  image  by  using  identical  pixel  coordinates.  Moving  target  signature  may 
exhibit  a  certain  amount  of  cross-range  coordinate  shifting,  defocusing,  and  smearing  in  the  SAR  image 
integrated  over  the  CPI,  but  otherwise  differ  from  stationary  clutter  by  virtue  of  the  phase  progression 
between  channels.  In  this  case,  SAR  images  from  multiple  channels  are  processed  to  calculate  the  AMF 
statistic  similar  to  Equation  (4. 1 ) 


samf  =  max 


"  vwR“'z 

2V 

vwR-'v 
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/ 

vn  =  un  exp 


jin 


nd  v, 
A  y 


a  J 


(4.2) 


Because  the  target  velocity  parameter  v,  is  unknown,  the  AMF  test  statistic  is  maximized  across  all 
possible  target  velocities.  Note  that.  Equation  (4.2)  is  the  same  as  Equation  (4.1)  except  for  a  functional 
change  in  v  by  invoking  the  matched  target  and  clutter  Doppler  condition,  vt/va=  cos#(sin$  -  sin0d). 

Although  the  single-cell  AMF  incorporating  “local”  array  calibration  results  in  excellent  clutter 
suppression,  this  approach  is  inapplicable  to  the  clutter-free  region  because  no  clutter  data  are  available 
for  array  calibration.  This  deficiency  motivates  an  alternative  “global”  calibration  technique  that  enables 
the  application  of  new  single-  and  multiple-cell  AMF  algorithms  to  both  the  endo-  and  the  exo-clutter 
regions.  These  AMF  processes  will  be  analyzed  in  detail  following  the  discussion  of  the  array  calibration 
procedure. 
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4.2  ARRAY  RESPONSE  CALIBRATION 


The  guiding  principle  of  correcting  the  array  distortion  is  to  derive  an  N  x  N  constant  matrix  C 
that  linearly  relates  the  actual  array  data  a((f>)  in  the  direction  (ft  to  the  ideal  array  response  v{(<f>)  as  a(tf>)  = 
C  vx((j>)  or  conversely,  Vi(<p)  =  C~ 1  a((p).  Because  numerous  clutter  cells  are  available  in  a  high-resolution 
range-Doppler  map,  one  may  accurately  estimate  C  via  the  standard  least-square  procedure  by  adopting  a 
measurement  data  model  X  =  CV  +  E.  The  data  matrix  X  =  [ a(<j) )),  ...  a(^)j),  ...]  consists  of  K  range- 
Doppler  clutter  data  vector  a(0/J,  k  =  1,  2,  ...  K,  V  =  ...  vx((j\),  ...]  is  a  similar  matrix  of  K  ideal 

Vandermonde  vectors,  vx(<f\),  k  =  1,  2,  ...  K,  and  the  Gaussian  matrix  E  models  the  measurement  errors. 
Assuming  an  aircraft  speed  va  and  a  depression  angle  #that  is  range-dependent,  the  azimuthal  angle  of 
a  range-Doppler  clutter  cell  is  related  to  its  Doppler/)  by/  =  (2v;/,^cos#sin^.  In  this  formulation,  the  C 
matrix  is  estimated  from  X as  follows: 


Array  calibration  model: 


a(<9)  =  Cv,(0)  or  v,(0)  =  C"'a(^) 
v,(0)  =  [l  ;  a  ;  a2  ...  /'']  ; 


;  C  is  0  dependent. 


f 


a  -  exp 


j27t—cos0sm(p 

A 


\ 


Measurement  data  model: 


X  =  CV  +  E  ; 


X  =  [a($)  a(<p2)  ...  a(^)  ...  a (<pK)]  ,  k  =  \,  2,  ...  K  ; 

V  =  [v,($)  v,(<z>2)  ...  v,((Z \)  ...  v,(^)]  ,  k  =  1,2,  ...  K 


LMS  estimateion  of  C:  C  =  XV  w  ( VVW ) 


(4.3) 


In  practice,  a  separate  calibration  matrix  C  is  estimated  for  each  range  gate  to  account  for  the 
elevation  angle  dependence  of  the  array  response.  This  technique  has  proved  to  be  very  effective  and 
satisfactory  in  correcting  the  array  response  distortion.  The  data  vector  a(<p)  in  Equation  (4.3)  also  may  be 
replaced  by  the  principal  eigenvector  of  the  local  clutter  covariance  R,  but  negligible  improvement  is 
expected  should  the  clutter-to-noise  measurement  be  high  enough. 

4.3  SINGLE-CELL  ADAPTIVE  MATCHED  FILTERING  WITH  “GLOBAL”  ARRAY 
CALIBRATION 

Once  completing  the  array  calibration  procedure,  the  single-cell  post-Doppler  AMF  detection 
amplitude  zAMF  is  now  calculated  from  Equation  (4.1)  by  using  the  calibrated  array  signal  v  =  Cv,  as 
derived  from  the  ideal  target  steering  vector  v, 
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v,(6>)  =  [l  ;  or  ; 


«2  ... 


j  ;  or  =  exp  y  2#— cos  6*  sin  $ 


(4.4) 


In  Equation  (4.4),  the  AMF  filter  output  is  specified  by  a  complex  amplitude  zAMF  instead  of  the 
power  iSamf  in  Equation  (4.1).  In  addition,  the  corresponding  minimum  variance  distortionless  response 
(MVDR)  filter  amplitude  zmvdr  is  specified  that  is  normalized  differently  as  an  unbiased  estimate  of  the 
signal  amplitude. 


4.4  MULTIPLE-CELL  SPACE-FAST-TIME  ADAPTIVE  MATCHED  FILTERING 


While  the  previous  single-cell  post-Doppler  AMF  or  MVDR  filter  is  capable  of  substantial  clutter 
suppression,  further  improvement  may  be  realized  by  passing  the  MVDR  output  through  a  Wiener  noise 
cancellation  filter.  The  Wiener  filter  (WF),  structured  as  a  sidelobe  canceller,  correlates,  estimates,  and 
partially  cancels  the  residual  interference  imbedded  in  the  MVDR  output  zMVdr  with  an  auxiliary  data  set 
y  that  assembles  range-delayed  copies  of  the  original  multichannel  range-Doppler  data  vector  x. 
Explicitly,  for  each  array  data  vector  x  at  a  range-Doppler  cell,  a  corresponding  y  data  vector  is  created  by 
concatenating  several  ^-vectors  of  adjacent  range  gates.  For  example,  given  M  range-delayed  taps 
centered  at  the  range-Doppler  cell  rd,  Mx  vectors  of  dimension  N  may  be  stacked  as  yvA  =  vec([*r.M/2  d ... 
*r-i,d  *r+i,d  VM/2,d])i  typically,  M  may  be  fixed  at  4  to  6  in  practical  data  processing.  Using  the 
auxiliary  data  y  of  dimension  MN,  the  detection  statistic  zWF  out  of  the  cascaded  MVDR  and  Wiener 
filters  is  derived  from  the  data  vector  z  of  the  test  target  cell  as 


v"R~'z 


»  ZWF  ~  ZMVDR  lx;yi'-yy-> 


R.-vR;y  ’ 


v  R  v 

R=vXxx"  ;  R,4iz“y"  ;  Rw=TZyyW  • 


K 

K I 

*•=]  k=\  ^  k= 1 

Auxiliary  data  vector  yrd  at  the  range-Doppler  cell  rd: 
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(4.5) 


yrd  =  vec 


d>~ 


Vi'd 


Xr-l,dXr+\,d 


The  procedure  used  for  estimating  the  covariance  matrices,  R:y  and  Ryy,  is  identical  to  the  previous 
scheme  used  for  estimating  the  covariance  R  of  x.  Note  that  the  adaptive  processing  DoF  is  now  equal  to 
the  dimension  of  they  vector  as  the  product  of  the  array  size  A  and  the  number  of  range-delay  taps  M. 

The  cascaded  MVDR  and  Wiener  filter  structure,  defined  by  Equation  (4.5),  actually  is  actually  a 
suboptimal  rendition  of  one  optimal  process.  The  optimal  process,  referred  here  as  multiple-cell  adaptive 
matched  filter,  swaps  the  order  of  execution  of  the  MVDR  and  the  Wiener  noise  cancellation  filters.  This 
novel  algorithm  is  applicable  to  any  situation  whenever  a  set  of  auxiliary  interference  plus  noise  data  y  are 
available  such  that  y  is  free  of  the  target  signal  and  is  correlated  with  the  primary  interference  process  x. 
In  the  current  application,  y  is  created  by  stacking  multiple  copies  of  the  range-delayed  array  data  vector* 
adjacent  to  the  test  target  range  cell  in  the  range-Doppler  domain  as  previously  explained.  Thus,  adaptive 
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processing  is  applied  to  the  two-dimensional  data  vectors  sampled  in  space  by  different  array  elements 
and  in  fast  time  by  different  range  gates. 

Assuming  that  the  auxiliary  data  y  may  suitably  be  modeled  as  a  complex  Gaussian  vector 
process,  one  may  derive  the  structure  of  the  optimal  multiple-cell  adaptive  matched  filter  and  its  detection 
statistic  by  performing  the  standard  hypothesis  test  on  the  composite  data  vector  zc  =  [z  ;  y\  to  determine 
the  presence  or  absence  of  a  target  in  the  test  data  vector  z  or  equivalently  zc.  Under  the  target  hypothesis 
Hi,  z  =  sv  +  x  is  the  superposition  of  both  a  target  signal  ,vv  and  the  interference  vector  x,  and  it  follows 
that  zc  =  (sv  +  x  ;  y].  Whereas  under  the  null  hypothesis  H0,z=x  consists  of  only  the  interference  x  and  zc 
-  pc ;  y].  Based  on  these  modeling  assumptions,  the  generalized  likelihood  ratio  test  computes  the  output 
detection  statistic  zm.-m  or  zAMr  of  the  multiple-cell  MVDR  or  AMF  algorithm  as  follows. 


Two-step  filtering  process: 

Wiener  noise  cancellation:  z„  =  z  -  R  R“'y  ; 

u  xy  yy  •' 


v"r:!z„ 


MVDR  filtering:  zMVDn  =  u_u,  u  ; 

v  R ...  v 


R  z 

or  AMF  zAMF  =  ■ - 


V^R 


LIU V  J 


R  is  the  covariance  of  residual  interference  u  after  Wiener  fitlering: 


u  =  x-RxyR;Jy  ;  Rull  =  Rxx  -  R^R^R^ 

K 


r,4Exx"  ;  Rxy=-^Xxyw  ;  Rw  =-4^  . 

A  *=1  A  k= I  A  *=l 


(4.6) 


The  multiple-cell  AMF  or  MVDR  processing  as  defined  by  Equation  (4.6)  uses  an  extended  array 
target  steering  vector  [v  ;  0]  of  higher  dimension,  whose  subspace  projection  v  is  identical  to  that  used  in 
the  single-cell  AMF  or  MVDR  filter.  Although  the  processing  steps  of  the  current  algorithm  are  well 
defined,  the  assumption  of  the  absence  of  any  target  signal  in  the  auxiliary  training  data  y  is  critical  and 
must  be  satisfied  to  achieve  high  interference  cancellation.  The  computational  load  of  the  multiple-cell 
algorithm  is  expected  to  be  higher  than  its  single-cell  counterpart  due  to  increasing  adaptive  processing 
DoFs,  now  equal  to  the  dimension  of  the  composite  data  vector  [x ;  y]  instead  of  the  x  data  vector. 

4.5  MULTIPLE-PIXEL  ADAPTIVE  SAR  IMAGE  PROCESSING 

A  similar  multiple-cell  noise  whitening  process  also  may  achieve  high  clutter  suppression 
performance.  This  scheme  is  motivated  by  the  realization  that  the  signature  of  a  moving  target  may  be 
smeared  after  integration  over  a  long  CPI.  The  signal  smearing,  often  neglected  in  the  preceding  single¬ 
cell  AMF  processing,  spreads  target  information  across  multiple  radar  range-Doppler  cells  or  SAR  range 
and  cross-range  pixels  in  each  channel.  For  this  reason,  adaptive  algorithm  performance  can  be  improved 
by  processing  the  extended  image  data  vectors  ze  constructed  by  concatenating  multiple  data  vectors  z  of 
adjacent  cells.  The  current  SAR  image  whitening  filter,  carried  out  in  the  extended,  multiple-pixel  vector 
space,  does  not  employ  a  target  model  but  incorporate  a  “noise  subtraction”  process: 

O  _7Wp-U 

°  whitening  ~  G-  ^  ^  e  ’ 

|(R)(  for  i  and  j  form  the  same  receiving  array  element  ;  (4.7) 

J  [  0  for  i  and  j  form  different  receiving  array  elements  . 
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Conceptually,  the  output  of  the  multiple-pixel  (or  multiple-cell)  whitener  consists  of  contributions 
from  (a)  same-pixel,  cross-element  terms,  (b)  cross-pixel,  cross-element  terms,  and  (c)  cross-pixel,  same- 
element  terms.  The  first  term  represents  information  produced  by  single-pixel  whitening.  The  second 
term  supplies  useful  new  information  about  moving  targets  not  available  with  single-pixel  whitening. 
The  third  pixel-space  whitening  term  does  not  carry  information  about  moving  targets  but  tends  to 
accentuate  the  noise  level  after  the  sample  covariance  matrix  inversion.  One  possible  remedy  is  to 
perform  noise  subtraction  by  employing  a  noise  covariance  R^that  is  constructed  from  R  by  zeroing  all  its 
cross-element  entries.  Consequently,  the  subtraction  in  Equation  (4.7)  should  cancel  the  purely  pixel- 
space  whitening  common  to  both  R  and  RN  while  leaving  the  cross-element  whitening  terms  intact. 
Mathematically,  one  may  regard  Equation  (4.7)  as  the  outcome  of  a  maximum  likelihood  ratio  test 
designed  to  test  the  presence  of  a  random  Gaussian  target  signal  vector  that  is  characterized  by  a 
covariance  matrix  R,  against  the  target-absent  data  that  are  instead  characterized  by  a  noise  covariance  RN. 

4.6  SAMPLE  RADAR  MEASUREMENT  RESULT 

One  sample  result  of  the  multiple-cell  AMF  incorporating  four  range-delayed  taps  is  depicted  in 
Figure  4.1  to  highlight  its  excellent  performance  in  clutter  cancellation.  The  data  used  in  this  example 
were  gathered  by  the  Tuxedo  radar  of  the  Lockheed  Martin  Corporation,  M&DS  Reconnaissance  System 
in  1999  under  the  DARPA  Discoverer  II  Program  [4],  The  Tuxedo  radar  antenna  has  three  horizontally 
polarized  subarray  apertures;  each  subarray  spans  0.61  m  in  length.  One  subarray  is  stacked  on  top  of  the 
other  two  adjacent  subarrays  so  that  the  total  azimuthal  aperture  length  is  1 .22  m.  Such  a  structure  yields 
the  azimuthal  beamwidth  of  3.6  deg  and  the  elevation  beamwidth  of  9.1  deg.  The  antenna  gain  is 
26.7  dBi.  The  particular  data  set  under  discussion  was  gathered  in  Pass  1916  by  the  Tuxedo  radar  in  its 
spotlight  mode.  The  radar  was  operated  at  the  center  frequency  of  9600  MHz  with  a  300  MHz  bandwidth 
and  the  radar  pulse  PRF  was  1363  Hz.  However,  the  data  bandwidth  was  reduced  to  180  MHz  for  this 
adaptive  processing  exercise.  In  addition,  the  data  pre-processing  procedures  also  included  channel 
equalization  over  the  operational  bandwidth  and  rudimentary  motion  compensation  to  correct  antenna 
phase  relative  to  a  fixed  aim  point  on  the  ground.  During  the  data  collection,  the  radar  flew  at  the 
approximate  altitude  of  5982  m  at  the  speed  of  136  m/s.  Five  ground  moving  vehicles  were  deployed  for 
measurement,  including  a  HMMWV,  a  2.5-ton  truck,  a  5-ton  truck,  a  large  fuel  truck,  and  a  PLS. 

Figure  15(a)  shows  the  range-Doppler  map  of  one  channel  after  coherent  integration  of  0.75  s. 
Prior  to  Doppler  filtering,  clutter  amplitude  has  been  limited  at  the  maximum  of  30  dB  above  noise 
without  alternation  of  its  phase.  The  range-Doppler  map  of  the  clutter  scene  generally  reflects  the  test 
ground  that  has  many  flat  surface,  roads,  and  fixed  structures  such  as  buildings.  The  map  indicates  the 
average  element  clutter-to-noise  ratio  within  the  mainbeam  region  is  below  20  dB.  However,  the  discrete 
clutter  strength  is  much  stronger  and  may  exceed  40  dB.  Figure  15(b)  plots  the  output  of  the  multiple-cell 
MVDR  process  from  Equation  (4.6)  that  depicts  the  signal-to-residual-clutter-plus-noise  ratio  in  the  same 
dynamic  range  and  over  the  same  target  area  as  that  in  the  Figure  15(a).  The  ability  of  this  adaptive 
process  to  cancel  ground  clutter  and  to  improve  the  test  target  signatures  is  apparent.  The  binary 
detection  map  of  the  Figure  15(c)  is  obtained  by  applying  a  sufficiently  high  detection  threshold  to  the 
output  MVDR  statistics  such  that  no  false  alarm  appears  in  the  data  domain.  In  this  example,  four  out  of 
five  moving  vehicles  are  detected.  A  range-cut  of  the  data  result  cross  the  slowest  target  at  the  range  of 
approximately  21.7  km  is  illustrated  in  the  Figure  15(d).  Note  that  the  slowest  target  being  detected 
moves  at  the  approximate  speed  of  0.5  m/s. 
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Figure  15.  Sample  output  of  the  multiple-cell  adaptive  processing  derived  from  the  Tuxedo  data  set  No.  1916 
over  a  test  area. 
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5.  SUMMARY 


A  unified  SAR  and  GMTI  radar  concept,  USASTAR  has  been  described  in  this  report.  USA 
STAR  provides  a  standoff,  rapid  wide  area  search  radar  concept  with  capability  to  detect  small  targets 
that  are  stationary  or  moving  at  arbitrary  speed  and  direction  simultaneously.  The  novel  concept  is  based 
on  large  baseline  multichannel  SAR  adaptive  clutter  cancellation,  while  minimizing  target  signal  losses. 
Target  signal  losses  are  minimized  by  operating  with  wideband  waveforms  and  antenna  baselines  that  are 
so  long  that  targets  move  two  or  more  resolution  cells  during  the  time  the  platform  flies  the  length  of  the 
antenna  baseline.  The  required  antenna  baseline  can  be  too  long  to  be  satisfied  with  filled  apertures,  but 
the  USASTAR  technique  would  work  with  very  sparse  arrays  such  as  two  separate  antenna  arrays  carried 
by  two  separate  platforms  flying  in  tandem,  as  long  as  phase  coherency  between  all  the  antenna  phase 
centers  is  assured  by  accurate  antenna  position  measurements. 

Theoretical  frameworks  formulating  the  USASTAR  concept  and  its  performance  are  provided  in 
this  report.  Performance  predictions  are  presented  in  terms  of  predicted  target  signal  loss,  or  SINR  loss, 
with  UHF  and  X-band  USASTAR  parameters.  Also  presented  are  algorithms  for  USASTAR  clutter 
cancellation  and  target  detection.  These  algorithms  were  successfully  tested  against  Tuxedo  multichannel 
X-band  data  on  moving  ground  targets.  Vehicle  targets  in  the  open  were  detected  at  a  radial  speed  as 
slow  as  0.5  m/sec.  Example  Tuxedo  X-band  results  are  presented. 
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